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We calculate the full transport counting statistics in a three-terminal tunnel device with one 
superconducting source and two normal-metal or ferromagnet drains. We obtain the transport 
probability distribution from direct Andreev reflection, crossed Andreev reflection, and electron 
transfer which reveals how these processes' statistics are determined by the device conductances. 
The cross-correlation noise is a result of competing contributions from crossed Andreev reflection 
and electron transfer, as well as antibunching due to the Pauli exclusion principle. For spin-active 
tunnel barriers that spin polarize the electron flow, crossed Andreev reflection and electron transfer 
statistics exhibit different dependencies on the magnetization configuration, and can be controlled 
by relative magnetization directions and voltage bias. 
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I. INTRODUCTION 

In crossed Andreev reflection (CA), a Cooper pair in 
a superconductor (S) is converted into an electron-hole 
quasiparticlc pair in normal-metal terminals (N„) or vice 
versa.i^ This process has potential applications in quan- 
tum information processing since it induces spatially sep- 
arated entangled electron-hole pairs. The performance of 
cntanglcrs utilizing this effect is diminished by parasitic 
contributions from electron transfers between the N n ter- 
minals. This process will be referred to as electron trans- 
fer (ET), but is also often denoted electron cotunncling 
or elastic cotunncling. 

It has been suggested that the noise properties 
of crossed Andreev reflection can be used to distin- 
guish it from electron transfer between the normal- 
metal terminals^ Mesoscopic transport noise also re- 
veals information about charge carriers which is in- 
accessible through average current measurements. In 
a system with several drain terminals, i.e. cur- 
rent beam splitters, noise measurements show signa- 
tures of correlations between the current flow in sep- 
arated terminals. The cross-correlation noise has at- 
tained interest recently since it can be utilized to 
study entanglement^^^^^^i 3 .^^^ and corre- 
lated transport. 

Crossed Andreev reflection in superconductor-normal- 
metal systems has been experimentally studied in Refs. 
[l7lll8lfl9l . In Ref. [l?] the nonlocal voltage was mea- 
sured in a multilayer Al/Nb structure with tunnel con- 
tacts between the superconducting Nb and the normal- 
metal Al layers. Current was injected through one of 
the normal-metal-superconductor contacts, and a non- 
local voltage measured between the superconductor and 
the other normal-metal. At injection bias voltage below 
the Thouless energy Et\i = HD/d 2 associated with the 
separation d between the normal- metals, positive non- 
local voltage was measured and this was interpreted as 
the result of dominating ET. For voltages eV above -&rh 
the nonlocal voltage changed sign, which was interpreted 



as a consequence of dominating CA. Subsequently, mea- 
surements reported in Ref. [l8| on nonlocal voltages in 
Au probes connected to a wire of superconducting Al by 
transparent interfaces indicated that the ET contribution 
is larger than the contribution from CA. 

The competition between CA and ET determines the 
sign of the nonlocal voltage and has been studied theoret- 
ically using various approache a 20 ' 21 i 22 ' 23 i 24 ' 25 i 26 ' 27 ' 28 in- 
cluding the circuit theory of mesoscopic superconductiv- 
ity utilized in this papcri 29 i 30 We will consider the linear 
response nonlocal conductance G n \. In supcrconductor- 
normal-mctal hybrid devices where transport in one 
normal-metal terminal Ni is measured in response to 
an applied voltage in another normal metal-terminal N2, 
this quantity is defined by 



dv 2 h — ~Gni — —{Get — Gqa), 



(1) 



where we have introduced conductances associated with 
the charge transfer processes introduced above, Gqa 
for crossed Andreev reflection, and Get for electron 
transfer^ The sign of the nonlocal conductance is de- 
termined by the competition between ET and CA. The- 
oretical calculations based on second order perturbation 
theory in the tunneling Hamiltonian formalism predicted 
that the nonlocal conductance resulting from CA reflec- 
tion is equivalent in magnitude to the contribution from 
ETi^i Thus the induced voltage in Ni in response to the 
bias on N2 should vanish since CA and ET give currents 
with opposite si gn, in contrast to the measurements re- 
ported in Refs. [l7l ll8|l9l. The tunneling limit was also 
considered in Refs. 26,29,31, and it has been found that 
the nonlocal conductance is in fact of fourth order in the 
tunneling and favours ET. 

Experimental investigations of crossed Andreev re- 
flection in superconductor-ferromagnet (S-F) structures 
have been reported in Refs. I32II33I . The measurements in 
Ref. were modeled using the theory of Ref. I20I . 

Experimental studies of the CA and ET noise proper- 
ties can be used to determine the relative contributions 
of these processes to the nonlocal conductance. It was 



2 



shown theoretically in that CA contributes positively to 
the noise cross correlations, whereas ET gives a negative 
contribution.— Calculations of higher order noise correla- 
tors or the noise dependence on spin-polarizing interfaces 
can reveal further information about the CA and ET pro- 
cesses. 

We will consider the full counting statistics (FCS) 
which encompasses all statistical moments of the current 
flowi2ii2£i2£ The noise properties of ET and CA reflec- 
tion thus obtained can be used to study the competition 
between these processes and reveal information that is 
not accessible in the mean currents. Our calculation also 
determines the contribution to the noise coming from the 
fermion statistics (Pauli exclusion principle). Moreover, 
the charge transfer probability distribution provided by 
FCS reveals information about the probability of elemen- 
tary processes in the circuit^ 

In this paper we calculate the FCS of multitcrmi- 
nal superconductor-normal metal and superconductor- 
ferromagnet proximity structures, and study the cur- 
rents, noise and cross correlations associated with the 
various transport processes. We obtain the probabil- 
ity distribution for transport at one normal-metal drain, 
and show that the probability associated with ET is 
larger than the probability associated with CA. For spin- 
active interfaces we show how spin filtering can be uti- 
lized to control the relative magnitude of the CA and 
ET contributions to the transport. Finally, we consider 
crossed Andreev reflection for spin triplet superconduc- 
tors. Thi s pape r goes beyond our previous publications 
Refs. Hiiflil in that we consider different bias volt- 
ages on the normal- metal/ferromagnetic drains and dis- 
cuss the effect of triplet superconductivity. 

The paper is organized in the following way: In Sec. |TT] 
we describe the electronic circuit and outline the formal- 
ism utilized to calculate the cumulant generating function 
of the probability distributions. In Sec. IIIII wc discuss the 
results in the case of normal metals, and in Sec. IIVI wc 
consider the spin-active connectors. Finally, our conclu- 
sions are given in Sec. [Vj 



II. MODEL 

The systems we have in mind can be represented by 
the circuit theory diagram (see Sec. Ill A|) shown in Fig. 
[TJ A superconducting source terminal (S) and normal 
metal drain terminals (N n ) are connected by tunnel bar- 
riers with conductance g n to a common scattering re- 
gion which is modeled as a chaotic cavity (c). The as- 
sumptions on the cavity is that the Green function is 
isotropic due to diffusion or chaotic scattering at the in- 
terfaces, and that charging effects and dephasing can be 
disregarded. The tunnel barriers can be spin-active with 
spin polarization <?MRn/ffn- We consider elastic trans- 
port at zero temperature. The superconducting terminal 
is grounded, and biases V n are applied to the normal 
terminals. We assume that V n <C Aq, where Aq is the 




FIG. 1: Circuit theory representation of the considered beam- 
splitter where supercurrent flows from a superconducting 
reservoir (S) into normal metal drains (N„). Tunnel barriers 
between cavity (c) and the drains can be spin-active, and are 
characterized by the conductances g n and spin polarizations 

gap of the superconducting terminal. In addition to the 
ET and CA processes described above, there can also be 
direct Andreev (DA) reflection between the superconduc- 
tor and one normal-metal terminal, where both particles 
of the Andreev reflected pair are transferred into N n . 
Semiclassical probability arguments show that the sub- 
gap charge current in the connector between Ni and c in 
the three terminal network in Fig. [1] has the following 
struct ure^^ 

h(E)= G CA (V 1 +V 2 )-G ET (V 2 -V 1 ) + 2G BA1 V 1 , (2) 

where we have introduced the conductance Gdai associ- 
ated with direct Andreev reflection between terminal N\ 
and S. Eq. |(2J) leads to the definition of the nonlocal con- 
ductance in ([1]) which shows that when V 2 > V\ , ET and 
CA give competing negative and positive contributions 
respectively to the current. The conductances in will 
be determined in the calculation below. 



A. Circuit theory 

The circuit theory of mesoscopic transport is reviewed 
in Ref. [jiH and is a suitable formalism to study proximity 
effects in superconducting nanostructures. The theory is 
developed from a discretization of the quasiclassical the- 
ory of superconductivity,— in combination with a theory 
of boundary conditions based on scattering theory. 

The circuit theory is formulated in terms of the qua- 
siclassical Green functions of the terminals and nodes in 
the system. Nodes can represent small islands or lat- 
tice points of diffusive parts of the system. Under the 
assumptions described above, the Green function of the 
spin-singlet S terminal in Fig. [I] is Gs = t\ where Tfe is a 
Pauli matrix in Nambu space. The Green functions G n 
of normal-metal terminals N„ are given by 

G (E) = I fsfs + + < eVn ' (S) 

c[ ' }f 3 f 3 + sgn( J B)f 3 (f 1 + l f 2 ) \E\>eV n , {) 
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where ffe arc Pauli matrices in Kcldysh space. 

Matrix currents / describe the flow of charge, energy 
and coherence between terminals and nodes through con- 
nectors, and conservation of these currents are imposed 
at each node. This generalized Kirchhoff law determines 
the Green function of the cavity G c . The balance of ma- 
trix currents flowing between each terminal n = 1, 2, S 
and the cavity, including the effect of superconducting 
pairing in c, can be written 



r c 2 



u Q V c 



A c fi, G c 



= 0. 



(4) 



Here, A c is the superconducting order parameter on the 
cavity, vo is the density of states and V c the volume of the 
cavity 30 i 39 The second term on the left hand side of the 
equation above induces electron-hole pairing. Since the 
pairing term has the same structure as the coupling to the 
superconducting terminal (see ©), it gives quantitative 
effects that are captured by renormalizing e 2 t/oV c A c /fi + 
.9s ~~ * .9s and will be omitted in the following. Here, gs is 
the superconductor tunnel conductance. 

Spin dependent transmission and reflection is de- 
scribed by tunneling amplitudes and r£ a for elec- 
trons with spin a incident from the cavity side on the 
interface between the cavity and terminal n in channel 
k. The matrix current /„ through such spin-active inter- 
faces is^ 



In — ~ [G„ , G c ] 



.9M Rn 



[{m„ • er r 3 ,G„} , G c ] . (5) 



Here, g„ = <?q ^ fc a |tj£ J 2 is the tunnel conductance 
where c/q = e 2 /h is the conductance quantum, gMRn = 
•9Q Sfcd^fe 1 1 2 — l^fe J 2 ) * s ^ e con ductance polarization, & 
is the vector of Pauli matrices in spin space, and the unit 
vector m„ points in the direction of the magnetization 
of the spin polarizing contact. In we have neglected 
an additional term related to spin dependent phase shifts 
from reflection at the interface which can be suppressed 
by a thin, non-magnetic oxide layer. 41 i 42 The effects of 
spin filtering contained in the polarization jMRn, which 
can be obtained experimentally using ferromagnetic ter- 
minals, will be studied in Sec. HVl 

In systems where all the connectors are tunnel barriers 
described by the matrix current ([5]) , it is possible to solve 
P| analytically and obtain the cavity Green function in 
terms of the terminal Green functions and the tunneling 
parameters. To this end, we note that it is possible to 
write ((4]) as [M, G c ] = 0. Employing the normalization 
condition G 2 = 1, the solution can be expressed in terms 
of the matrix M asii 



G c 



(6) 



This result facilitates calculation of the cumulant gener- 
ating function of the charge transfer probability distribu- 
tion in tunnel barrier multitcrminal circuits. 



B. Full counting statistics 

Full counting statistics is a useful tool to compute 
currents and noise in a multiterminal structure^ 3 - and 
also provides the higher statistical moments that may 
become experimentally accessible in these systems. Ad- 
ditionally, one can obtain information about the elemen- 
tary charge transport processes by studying the proba- 
bility distributions^ The cumulant generating function 
(CGF) S({\ n }) of the probability distribution is directly 
accessible by the Green function method, and is defined 
by 



-5({x„}) 



:^P({iY„};t )e- i ^ 



XnN n 



(7) 

Here, P({2V„};io) is the probability to transfer 
Ni, N2, ■■■jNn electrons into terminal Ni,N2, ...,N„ in 
time to, and M is the total number of terminals in the 
circuit. The CGF is a function of the set of count- 
ing fields {xn} which are embedded in the Green func- 
tion at each terminal by the transformation G n — > 
e i Xn fK/2Q ne -ix n fK/2 w here f K = f 3 fi. The CGF will 
be determined by the following relation,— 



ie dS({xn}) 
to dxn 



dEI n ({ X n}), 



(8) 



where I n {{Xn}) is th e particle (counting) current through 
connector n in presence of the counting fields. Our task 
is now to integrate this equation and obtain the CGF 
S({xn})- Using the general solution to the matrix cur- 
rent conservation ([6]), it was found in Ref. [Ill, that this 
is possible by rewriting the counting current in terms of 
a derivative of M with respect to the counting fields. 
Explicit derivation shows that 

4^ x „Tr{^}. 

(9) 

This result is valid also in the presence of spin-active 
contacts ©. Combining © with © yields S({x n }) 
straightforwardly. 

Practical calculations of CGFs are performed by diag- 
onalizing the matrix M, which allows us to express the 
CGF in terms of the eigenvalues of M, 



M{Xn}) = ^Tr{fic/„({x4)} 



S 



to_ 

"4e 2 



dE 



E 

fe 



(10) 



In this equation, is the set of eigenvalues of M. 

We can obtain the cumulants of the transport proba- 
bility distribution by successive derivatives of the CGFi2£ 
Specifically, we obtain the mean current from 



^ = _te dS({ X }) 

to dXn 



(11) 



{x=o} 
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The current noise power is given by 
^ d 2 S({ x }) 



2 G - 



to dXmdXr 



(12) 



{x=o} 



where in the multiterminal structure, the autocorrelation 
noise at terminal n is given by C„ iri . When m n, (TH 
gives the noise cross-correlations. 



III. NORMAL METAL DRAINS 

In this section, we will consider the FCS of the super- 
conducting beamsplitter in Fig. [1] when the connectors 
are not spin polarizing, and generalize previous works by 
taking into account a difference in drain terminal voltages 
Vl ^ V 2 . 

In the regime E < eV% , eV 2 the only contribution to the 
nonlocal conductance comes from CA since we consider 
zero temperature. The resulting CGF was studied in Ref. 
UTI where it was assumed that V\ = V%. In the general 



J 



case Vi ^ V 2 , the total CGF S following from (10} has 
one contribution from the energy range E < eV\,eV2, 
and if the voltages are different, another contribution in 
the energy range eV\ < E < eV 2 (we assume V 2 > Vi), 




sM+SbiVi-Vi). 



(13) 



There is no contribution to transport at \E\ > eV 2 - Here, 
we have defined two separate contributions to the CGF 
that govern transport in the regime E < eVi,eV 2 (iS a ) 
where only Andreev reflections (CA and DA) can occur, 
and the regime eV\ < E < eV 2 (5b) where in addition to 
Andreev reflections, ET can take place. The contribution 
S a was calculated in Ref. [ll|, see (|14a[) . where we have 
defined gs = [g$ + g 2 ] 1 ^ 2 and g ~ gi + .92- The counting 
factors e^xs-ixm-ixn describe processes where two par- 
ticles are transferred from S, and one particle is counted 
at terminal N TO and at terminal N„ (m, n = 1, 2). 



5 a 



s h = - 



1 



tpVi 



\ 



gl + g* + 4g 2 ^ g m g n {^x^x m -^ n - 1), 



V2e 



(14a) 



gl + 2g l92 (e^-^ -1 




4 1 a„2S^ g n g 2 ( e 2i Xs -ixn-iX2 - 1) + 4g 1 g 2 g 2 (e i ^- i ^ - 1) 



(14b) 
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FIG. 2: Transport processes in the three terminal device when eVi = 0: (a) Crossed Andreev reflection: A Cooper pair from S 
is converted into an electron-hole pair in c by Andreev reflection, and the electron with energy +E is transferred into Ni, and 
the hole with energy — E is transferred into N2. Tunnel barriers between the reservoirs may be spin-active and are described by 
magnetization vectors m that in this paper are considered collinear. (b) Electron transfer: A particle from Ni tunnels through 
the cavity c into N2. The density of states in the cavity c is suppressed due to the proximity effect from the superconducting 
terminal. 



In (|14b[) we show the calculated 5b which has con- tributions from electron transfer. The counting factor 
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c ?Xl ~* X2 describes events where an electron is transferred 
from Ni to N 2 . Compared to 5 a , we see that DA events 
between S and Ni that would be described by counting 
factors e 2iXS_2lXl no longer occur. This can be under- 
stood from the electron-hole-nature of Cooper pairs, see 
Fig. [2l Two quasiparticles, with energy +E for the elec- 
tron and — E for the hole constitute the Andreev reflected 
quasiparticle pairs in c. In the energy range considered 
here, eVi < E < eV 2 , the states in Ni are occupied, pre- 
cluding DA reflection into Ni . A similar argument shows 
that in CA processes, the electron must be transferred 
into Ni and the hole into N 2 . 

The nonlocal conductance G n i = GW — Gca following 
from (|14b[) is in agreement with Ref. [29| where 



G 



ET 



5152 2 5 2 



9l 



2 [<? 2 + 5§] 3/2 



, G 



CA 



5i52 



2 [5 2 +5i] 3/2 ' 



(15) 



The nonlocal conductance is dominated by ET and is 
of order 0{g^)M^ When g/g s <C 1, the nonlocal con- 
ductance vanishes due to equal probability for ET and 
CA as we will explicitly show by inspection of the prob- 
ability distribution below. In the opposite limit that 
the coupling to normal terminals dominates, 5s/5 <C 1, 
the conductance for ET is to lowest order given by 
5i52(l - SN)/g. Here SN = (g s /g) 2 /2 is the lowest or- 
der correction to the low-energy density of states due to 
superconducting correlations. The CA conductance be- 
comes gig2SN/g in this limit. 

Let us now consider cross-correlations. In general, CA 
leads to a positive contribution and ET leads to a nega- 
tive contribution to the cross-correlation. An additional, 
negative contribution is induced by the Pauli exclusion 
principle. The cross-correlation between Ni and N 2 fol- 
lowing from (|13p is 



Ci, 2 = 2e(Vi + V 2 )G CA - 2e(V 2 - Vi)G ET (16a) 

_I^( Gl -G nl )(G 2 -G nl ) (16b) 
5s 

4e(V a -Vi) 



5s 



[Gca(G 2 + 2GdA2) - GdA2G„i] , (16c) 



where we have defined the local differential conductances 
G„ = d n I n and the conductance for direct Andreev re- 
flection into terminal n is GuAn = 5s5 2 /25e^ We now 
focus on the competition between CA and ET. When the 
two normal-metal terminals are at equal voltage V\ = V 2 
the contribution from ET and also the term in (|16c|) van- 
ishes and we are left with a positive contribution to the 
cross-correlations from CA due to the correlated parti- 
cle transfer into Ni and N 2 . An additional, negative 
contribution due to the Pauli principle in (|16b[) van- 
ishes in the limit of asymmetric conductances gs 3> g or 
5i(2) 3> 52(1) j 5s due to the noisy (Poissonian) statistics 
of the incoming supcrcurrent. A negative contribution 
from ET in (|16ap is induced when there is a voltage dif- 
ference between the normal-metal terminals due to the 



currents with opposite signs in Ni and N 2 resulting from 
this process. This demonstrates that it is possible to tune 
the sign of cross-correlations by the voltages V\ and V 2 .— 
The contribution to Ci, 2 in p6c[) is proportional to the 
voltage difference V\ — V 2 , and vanishes in the limit of 
asymmetric conductances gs ^ g or 5i( 2 ) S> 5 2 (i),5s- 

It is interesting to compare 5b with the corresponding 
CGF when S is in the normal state, 



Sh = 2e X 



(17) 



V(9i + 5s + 52) 2 + 4g 2 gie lx i-' x 2 + Ag 2 g S Q ixs ~ iX2 ■ 

Here we see a contribution due to transport between 
Ni and N 2 that is similar to the one outside the dou- 
ble square root in ()14b[) . Superconductivity leads to the 
double square root in (|14bp that takes into account the 
correlation of transport through c by Andreev reflections 
and ET. The complicated dependence on the counting 
fields in (|14bj) precludes a simple interpretation of 5b in 
terms of the probabilities of elementary charge transfer 
processes. However, when 5s 3> 5 or g\(2) 3> 52(1)) 5S we 
can expand the square roots in 5b and obtain the CGF 



5 b 



t (V 2 -Vl) 

25|e 



g ig2 (g 2 s + 2g 2 )e^-^ 
+ g 2 gi g 2 e 2 ^-^-^ 



(18) 



In this limit the CGF describes independent CA, ET, 
and DA Poisson processes. The prefactors determine the 
average number of charges transferred by each process in 
time to- 

To illustrate the physics described by 5b in the limit 
introduced above, let us examine the probability distri- 
bution obtained by the definition in ([7]). If we consider 
the current response in Ni to a voltage in N 2 , we can 
consider that Vi = and the only contribution to the 
total CGF 5 comes from 5b- The normalized probability 
distribution for the transport at terminal Ni following 
from ([TBI) then becomes 



P(N i; t Q ) =e 



E 



fe>|JVi| 
k-\-Ni even 



N 4 



Here we have defined the mean number of particles trans- 
ferred in time to, 

\h\to V 2 t 



51.92 



e g{l+g 2 /g 2 f 2 ' 



(20) 



Eq. (flT)]) describes a joint probability distribution for CA 
and ET processes with Poissonian statistics, and is con- 
strained such that the number of CA events described 
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by the weight g^/2g 2 subtracted by the number of ET 
events described by the weight ffg/2g 2 + 1, is Ni as re- 
quired. When gs/g 3> 1, the mean number of particles 
transferred vanishes according to (j2"0|) and the probabil- 
ity distribution (|19[) is symmetric around Ni = 0. This 
means that the average current vanishes, since the prob- 
abilities for ET and CA are equal. In general, the proba- 
bility distribution has its maximum for negative N±, i.e. 
ET is more probable than CA reflection. In Fig. [3] wc 
have plotted the probability distribution for differ- 
ent values of gs/g- For small ratios gs/g, ET dominates 
and the probability distribution is centered at a nega- 
tive value for iVi. As expected, we see that the center 
of the probability distribution (mean number of parti- 
cles transferred) is shifted from a negative value towards 
zero with increasing gs/g- The width of the distribution, 
described by the autocorrelation noise C\,i (see (|12|)). 
decreases with increasing gs/g- 



IV. SPIN-ACTIVE CONNECTORS 



Qualitatively, the effect of spin polarizing interfaces 
on the competition between CA and ET processes in S-F 
systems can be understood as follows^ The ET process 
is favoured when magnetizations of ferromagnetic leads 
are parallel since the same spin must traverse both the 
interfaces between c and the ferromagnets. On the other 
hand, CA reflection is favoured in an antiparallel config- 
uration since two particles with opposite spins must tra- 
verse the interfaces. This behaviour was experimentally 
observed in Refs. I32ll33l where ferromagnetic Fe probes 
were contacted to a superconducting Al wire. 
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The FCS of a beam splitter with spin-active contacts 
and V\ = Vi was considered in our previous paper Ref. 
l37l and constitutes iS a , see (|21a|) . In this case the only 
transport processes are DA and CA reflection, and we 
found that CA is enhanced in an antiparallel alignment 
of the magnetizations as expected. 



FIG. 3: (Color online) Probability distribution for transport 
of Ni electrons into terminal Ni, P(Nx;to) Eq. [fl9j). We 
show distributions for two different values of the parameter 
gs/g = 5.00 (red [gray] solid impulses), and 0.05 (blue [dark 
gray] dotted impulses). We have chosen the parameter a — 
gig2V2to/eg — 20, which gives the mean value Ni for small 
g s /g, see Eq. (g0j|. 



When V2 > V±, there can also be ET, and an addi- 
tional effect is spin accumulation at the node. With 
collinear magnetizations (sign of gyiRn describes mag- 
netization directions up (positive) or down (negative) 
along the z quantization axis) , we find that in the regime 
eVi < E < eV2, Sb = Sb+ + <Sb- where iSb CT is given 
below (plb|) . 
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t v 

V2e 



\ i 



4 ffMR^2 , , 4 5 2 



{gl + .9 2 ) + -rXi (5mffn ~ SMRmSMRn) (e«X-<Xm-ix» - 1), 



(21a) 



5 b 



2V2e I ff E 

1 - ^¥(9 2 s + 9 2 ) + A ~4 £(ff» - ^MR B )fe + <7<7mr2 ) (e 2 ^-^-^ - l) 
9s 9t. „ 



4(.g - 3mr) 5 



9i 



(.91 + 0-ffMRl)(52 + 0-ffMR2) (c* Xl 4X2 - 1) 



-,1/2' 



1/2 



(21b) 



Here we have redefined g% = (3s+3 2 +5mr) 1 ^ 2 an d intro- 
duced jmr = <7mri + <?mr2- The expression for St reduces 
to the result for nonpolarizing contacts, Eq. (|14b[) . in the 
limit that jMRn ~> 0. The two terms S^a correspond to 
the two possible directions of the spin(s) involved in the 
charge transfer processes. The spin-dependent conduc- 
tance for a spin up (down) is g n + (— )<?mr«- In ET, one 
spin must traverse the two spin-active interfaces, thus the 
counting factor for spin a is proportional to the weight 
(gi + <7<?mri ) (52 + c<?MR2)' The two spin channels are 
independent. The two opposite spins of an Andreev re- 
flected quasiparticle pair can be CA reflected into termi- 
nals with different polarizations <?mri and (?mr2 accord- 
ing to the prefactor (gi — ct5mri)(<?2 + CTC/MR2), an d each 
possibility for the directions of the two spins gives an 
independent contribution to the CGF 5b- 

In the limit that g s > g or g 1(2) > .92(1), ffs (<7MRn < 
g n by definition) we can expand the double square roots 
and perform the summation over Shu which yields 



5 b 



t Q (V2-Vi) 



x { [fll + (.9 - 3mr) 2 ] (gm + ffMRiffMR 2 )e lXl 4X2 
+ 9 2 s (.91.92 - .9MRi.9MR 2 )e 2lxs -* Xl - lX2 
+.9i(.92 2 -5 2 , R2 )e 2lxs - 2lX2 }. (22) 
The nonlocal conductance following from l]22[) is given by 

[gl + (g - !7mr) 21 



Get —(51.92 + .9mri5mr2) 



GcA =(.9lff2 — 3MR.1.9MR.2) 



2<?r 



(23a) 
(23b) 



This immediately demonstrates that ET is favoured 
in a parallel configuration of the magnetizations 
(<7mri3mr2 > 0) as the same spin in this case tun- 
nels through both interfaces. On the other hand, 
CA reflection is favoured by antiparallel magnetizations 
(<?mri5mr.2 < 0) since the opposite spins of a singlet tun- 
nel through different interfaces. These qualitative fea- 
tures are in agreement with Ref. I20I . 



The cross-correlation following from (|22l is 
Ci,a = 2e(y 2 + Vi)Gca - 2e(V 2 - Fi)G ET . 



(24) 



The sign of C12 can now be tuned by two experimental 
control parameters: The bias voltages through the pref- 
actors in (|24p . and the relative magnetization direction 
that determines the magnitudes of Get and Gqa- 

In the energy range V\ < E < V% we are in this setup 
measuring the energy of the quasiparticlcs involved in 
crossed Andreev reflection, see Fig. [2j Since the electron- 
like quasiparticle with energy +E must flow into Ni, 
and the —E hole-like quasiparticle must flow into N2, 
this means that the entanglement in the energy degree 
of freedom of Andreev reflected quasiparticle pairs has 
collapsed. 



A. Triplet superconductivity 

Superconducting correlations with triplet pairing sym- 
metry in the spin space can be induced by magnetic 
exchange fields in singlet superconductor heterostruc- 
tures. This effect has attained considerable interest, and 
has recently been experimentally demonstrated (see Refs. 
I44ll45ll46l and references within) . The different spin-space 
symmetry of the induced electron-hole correlations opens 
interesting experimental applications where e.g. super- 
conducting correlations can propagate through a strong 
ferromagnetic materiaL 47 i 48 We have studied the FCS 
when S is a source of spin triplet quasiparticle pairs, and 
in this subsection we summarize our results for collinear 
magnetizations when V\ =V-2- 

The CGF for the spin triplet Cooper pairs with 
S z I "J) = 0, where S z is the spin operator along the z- 
axis and is the spin part of the Cooper pair wave 
function, is identical to the CGF for conventional spin 
singlet superconductors. We showed in Ref. H3 that the 
CGF (|21ap reveals the entangled nature of the quasipar- 
ticle pairs. The S z \^>) = spin triplet states are also 
one of the maximally entangled Bell states which implies 
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that the magnetization dependence for CA is the same 
for singlet and triplet in the collinear case. This result 
can be shown also straightforwardly by computing the 
two-electron tunneling probability which is propor- 
tional to (V\ (51 + ffMRlCT 3 ) ® (92 + 5MR2^s) |*) + 1 *-> 2. 



The triplet states where S z = ±fi |\&) give rise to a 
different dependence on the magnetization configurations 
in the CGF since the quasiparticle pairs are not in spin 
entangled states, but rather in product states. We obtain 
the CGF 



t V 
2v / 2e\ 



4.9mr 



■J MR 

9i 



(fl§ +9 2 ) + ^r {91 + <J9mki) {g-2 + ag M R2) (e 2 ^-ix™-ixn _ l), 



9i 



(25) 



Compared to the singlet case (|21a[) . we see that the 
CA counting prefactor factorizes as a result of the non- 
entangled product state of the quasiparticle pairs. This 
magnetization dependence can be recovered also by cal- 
culating the two-particle tunneling probability p\ 2 as dis- 
cussed above. 



V. CONCLUSION 

We have calculated the full counting statistics of mul- 
titerminal tunnel-junction superconductor-normal-metal 
and superconductor-ferromagnet beam splitter devices, 
and studied the resulting currents and cross-correlation. 
The probability distribution for transport at a normal- 
metal drain contact demonstrates that the probability 
for electron transport between two normal-metal termi- 
nals is larger than the probability for crossed Andreev 
reflection. A finite voltage difference between the normal- 
metal contacts introduces competing contributions to the 



cross-correlations from electron transport between nor- 
mal terminals and crossed Andreev reflection. Finally, 
we have shown how spin-active contacts act as filters for 
spin, and calculated the cumulant generating function. 
The sign of the cross-correlation due to the competing 
contributions from electron transport between drain ter- 
minals and crossed Andreev reflection can in this case 
be determined by two external control parameters, i.e. 
bias voltages and the relative magnetization orientation. 
Finally, we make some remarks about the counting statis- 
tics in the case of spin triplet superconductors. 
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